Multivariate piecewise linear interpolation of a random field 
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Abstract 

We consider a multivariate piecewise linear interpolation of a continuous random field on a d-dimensional 
cube. The approximation performance is measured by the integrated mean square error. Multivariate 
piecewise linear interpolator is defined by N field observations on a locations grid (or design). We investigate 
the class of locally stationary random fields whose local behavior is like a fractional Brownian field in 
mean square sense and find the asymptotic approximation accuracy for a sequence of designs for large N. 
Moreover, for certain classes of continuous and continuously differentiable fields we provide the upper bound 
for the approximation accuracy in the uniform mean square norm. 
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Let a random field X(t), t € [0, l] d , with finite second moment be observed at finite number of points. Suppose 
further that the points are vertices of hyperrectangles generated by a grid in a unit hypercube. At any 
unsampled point we approximate the value of the field by a piecewise linear multivariate interpolator, which 
is a natural extension of a conventional one-dimensional piecewise linear interpolator. The approximation 
accuracy is measured by the integrated mean squared error. This paper aims m odelling random fields with 
given accuracy based on a finite number of observations. Following iBermanl (|1974l ). we extend the concept of 
local stationarity for random fields and focus on fields satisfying this condition. For quadratic mean (q.m.) 
continuous locally stationary random fields, we derive the exact asymptotic behavior of the approximation 
error. A method is proposed for determining the asymptotically optimal knot (sample points) distribution 
between the mesh dimensions. We also study optimality of knot allocation along coordinates of the sampling 
grid. Additionally, for q.m. continuous and continuously differentiable fields satisfying Holder type conditions, 
we determine asymptotical upper bounds for the approximation accuracy. 

The problem of r andom field approximation arises in many re search and applied areas, like Ga ussian ran 



dom fields modell ing (jAdler and Taylor]. 



1992 



2007 



Brouste et al 



Stein, 1999), sensor networks (I Zhang and Wicked . 120051 ) , and image processing (jPrattl . 120071 ) . The upper 



20071 ) , environmental an d geoscience s (jChristakoi 



bound for the app ro ximation error for isotro pic random fields satisfying Holder type conditions is given in 



Ritter et al. 



|l995|). 



Miiller-Gronbach 



(|l998l ) consider affine linear approximation methods and hyperbolic 
cross designs for fields with covariance function of tensor type. An optimal allocation of the obseryation s 
for Gauss ian random fields with product type kernel is investigated in lMiiller-Gronbach and Schwabd (|1996l ). 



Su 



(|1997l ) studies limit behavior of the piecewise constant estimator for random fields with a particular form 

) investigates exact asymptotics of stationary spatial process approx- 



Benhenni (2001 



of covariance function. 

imation based on an equidistant sampling. The approximation complexity and the curse of dimensionality 
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(2003); lAbramowicz and Seleznjevl (120111 ). iRitterl (120001 ) contains a very detailed survey of various random 



in, e.g., 


Selezniev 


( 


1996). Re- 


Seleznie\ 


( 


2000): 


Hiisler et al. 



process and field approxi mation pr o blems . For an extensi v e stud i es of approximatio n problems in deterministic 



setting, we refer to, e.g., iNikolskiil (jl97fil ): Ide Boor et al.1 (|2008l ): IKuo et all (120091 ) 



The paper is organized as follows. First we introduce a basic notation. In Section 2, we consider a piecewise 
multivariate linear approximation of continuous fields which local behavior is like a fractional Brownian field 
in mean square sense. We derive exact asymptotics and a formula for the optimal interdimensional knot 
distribution. In the second part of this section, we provide an asymptotical upper bound for the approximation 
accuracy for q.m. continuous and differentiable fields satisfying Holder type conditions. In Section 3, we present 
the results of numerical experiments, while Section 4 contains the proofs of the statements from Section 2. 

1.1 Basic notation 

Let X = X(t),t G £> := [0, l] d , be a random field defined on a probability space (O, J£~, P). Assume that for 
every t, the random variable X(t) lies in the normed linear space L 2 (Q) = L 2 (ft, & , P) of random variables 

1 /2 

with finite second moment and identified equivalent elements with respect to P. We set ||£|| := (E£ 2 ) 
for all £ G L 2 (£l) and consider the approximation based on the normed linear spaces of q.m. continuous and 
continuously differentiable random fields denoted by C{T>) and C l (T>), respectively. We define the norm for 
any X £ C(V) by setting 



XI 



\X(t)\\ p dt 



v 



i/p 



1 < p < CO, 



and || X := maxtg© || X(t) ||. For p = 2, we call the norm integrated mean squared norm and the corre- 
sponding measure of approximation accuracy the integrated mean squared error (IMSE). 

Now we introduce the classes of random fields used throughout this paper. For k < d, let 1 = (1%, 



l h i ' 



0, 



k, Ln 



> h) 
0, be the 



be a vector of positive integers such that Ylj=i h = d, and let L, := Y^j 
sequence of its cumulative sums. Then the vector 1 defines the I- decomposition of T> into V 1 x T> 2 x . . .V k , 
with the Zj-cube T>i = [0, lp , j = 1, . . . , k. For any s £ D, we denote the coordinates vector corresponding to 
the j'-th component of the decomposition by s J , i.e., 



s i = s i(l):=(s L . 1+1 ,...,s L .)eVi, j = 1, . 



k. 



For a vector a = (a±, . . . , «&), < ctj < 2, j = 1, . . . , k, and the decomposition vector 1 
define 



(h, ■ ■ ■ ,h), w e 



3=1 



with the Euclidean norms | |s J 1 1 , j = 1, . . . , k. 



j \\ j forallseP 



For a random field X G C([0, l] d ), we say that 

(i) X G CpQO, l] d , C) if for some a, 1, and a positive constant C, the random field X satisfies the Holder 
condition, i.e., 



X(t + s) -X(t)|| 2 < C||s| 



for all t,t + s G [0,l] d , 



(1) 
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(ii) X G #j*([0, l] d ,c(-)) if for some a, 1, and a vector function c(t) = (ci(t), . . . , Cfc(t)), t G [0, l] d , the random 
field X is locally stationary, i.e., 

||X(t+s)-X(t)|| 2 ^ uniformly in t G [0, l] d , (2) 

E-=ic fc (t)||^||^ 

with positive and continuous functions ci(-), . . . , Cfc(-). We assume additionally that for j = l,...,/c, the 
function Cj(-) is invariant with respect to coordinates permutation within the j-th component. 



For the classes C" and Bf, the withincomponent smoothness is defined by the vector a = (oti, . . . , a^). 
We denote the vector describing the smoothness for each coordinate by a* = (a*, . . . , a£), where a* = ay, 
i = Lj-i + 1, . . . ,Lj, j = 1, . . . , k. 

Example 1. Let m = (mi,...,mfc) be a decomposition vector of [0, l] m , and m = Y^]=\ m r Denote by 
-B/3,m(t), t G [0, l] m , (3 = (Pi, . . . ,/3fc), < /3j < 2, j = 1, . . . , k, an m-dimensional fractional Brownian field 
with covariance function r(t,s) = ^ (| |t| 1^ + | |s||jg — ||t — s\\p). Then Bp m has stationary increments, 

||^/3, m (t + s)-£? Am (t)|| 2 = ||s|l/3, t,t + sG [0,l] m , 

and therefore, Bp m G B^V,^-)) with local stationarity functions ci(t) = . . . = c^(t) = 1, t G [0, l] m . In 
particular, if k = 1, then Bp >m (t), t G [0, l] m , < /? < 2, m G N, is an m-dimensional fractal Brownian field 
with covariance function 

r(t,s) = i(||t|| /3 + ||s|| /3 -||t-s|| /3 ), t ,t + sG [0,l] m . (3) 



For X G ^([0, l} d ), we write X'^t), t G [0, l] d , to denote a q.m. partial derivative of X with respect to the 
j-th coordinate, and say that X G C 1 '"* ([0, l] d , C) if there exist a vector a* = (aj,...,a£) and a positive 
constant C such that each partial derivative Xj is Holder continuous with respect to the j-th coordinate, i.e., 
if for allt,t + sG [0, l] d , 

WX'jih, tj+Sj, . . . ,t d ) - X'jih, . . . ,tj, . . . , t d )\\ 2 < C\sj\ a i j = l,...,d. (4) 

Moreover, we say that X G Cj 1,a ([0, l] d , C) with a = ( ai ,...,a k ) if X G C 1 '"* ([0, V\ d , C) and for a given 
partition vector 1, ctj := ol* l +1 = . . . = a* L ., i = 1, . . . , k. 

Let X be sampled at N distinct design points T/v- We consider cross regular sequences of sampling designs 
T/v := {ti = (ti^j , . . . , t^i d ) ■ i = (i\, ■ ■ ■ , id), < < n* k , k = 1, . . . , d} defined by the one-dimensional grids 

f tj ' i i 

/ h*(v)dv = —, i = 0,l,...,n*, j = l,...,d, 
Jo n j 

where h*(s), s G [0,1], j = l,...,d, are positive and continuous density functions, say, withindimensional 
densities, and let 

h*(t) := {ht(h),...,h* d (t d )). 
The interdimensional knot distribution is determined by a vector function it : N — > ¥\ d : 

vr*(iV) :=(nl(N),...,n d (N)), 
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where hnijv_ s . 00 n*(N) = oo, j = 1, . . . , d, and the condition 

d 

H(n*(N) + l) = N 
3=1 

is satisfied. We suppress the argument N for the sampling grid sizes n*j = n*(N), j = 1, . . . , d, when doing 
so causes no confusion. Cross regular sequences are one of the possible extensions of the well known regular 



sequences introduced by ISacks and Ylvisakerl (| 19661 ) . The introduced classes of random fields have the same 
smoothness and local behavior for each coordinate of components generated by a decomposition vector 1. 
Therefore in the following, we use only approximation designs with the same within- and interdimensional 
knot distributions within the components. Formally, for the partition generated by a vector 1 = . . . , Ik), we 
consider cross regular designs Tv, defined by the functions h := (hi, ...,%) and it(N) := (ni(N), . . . , rik(N)), 
as follows: 

K(-) = h j(-)^ n i= n ji i = Lj-i + l,...,Lj, j = l,...,k. 

We call the functions hi(-), . . . , %(•) and ir(N) withincomponent densities and intercomponent knot distribu- 
tion, respectively. The corresponding property of a design T/v is denoted by: Tv is cRS(h,ir,Y). 

For a given cross regular sampling design, the hypercube T> is partitioned into M = Ilj=i n j disjoint 
hyperrectangles T>\, i = (ii, . . . , id), < ik < ra£ — 1, k = 1, . . . , d. Let 1^ = (1, . . . , 1) denote a (i-dimensional 
vector of ones. The hyperrectangle T>\ is determined by the vertex t; = (ti i 1} . . . ,tji d ) and the main diagonal 
*"i = t i+ i d - t i; i.e., 

D, : = |t : t = ti + ri*s,s G [0, , 

where denotes the coordinatewise multiplication, i.e., for x = (x\, . . . ,Xd) and y = (y\, . . . ,ya), x * y := 
{xiyi, . . .,x d y d ). 

For a random field X G C(T>), define a multivariate piecewise linear interpolator (MPLI) with knots T/v 

X N (t):=X N (X,T N )(t)=E v X{t i + r i *r 1 ), t G V h t = t; + r ; * s, 

where r\ = (rji, . . . , rjd) and rji,...,rid are auxiliary independent Bernoulli random variables with means 
si,...,Sd, respectively, i.e., r\j G Be(sj), j = l,...,d. Such defined interpolator is continuous and piece- 
wise linear along all coordinates. 

Example 2. Let d = 2, N = 4, V = [0, l] 2 . Then t = s, r = (1, 1), 

Xjv(t) = Vr,X(r]) = X(0, 0)(1 - h)(l - t 2 ) + X(l, 0)h(l - t 2 ) + X(0, 1)(1 - h)t 2 + X(l, l)ht 2 , 



and Xtv is a conventional bilinear interpolator (see, e.g.. lLancaster and Salkauskasl . 11986). 

We introduce some additional notation used throughout the paper. For sequences of real numbers u n 
and v n , we write u n < v n if lim n ^. 0O u n /v n < 1 and u n x v n if there exist positive constants c\,c 2 such that 
c\u n < v n < c 2 u n for n large enough. 

2 Results 

Let Bp tTn (t),t G W+, < (3 < 2, m G N, denote an m-dimensional fractional Brownian field with covariance 
function ([3]). For any u G M™, we denote 

bp, m (u) := / || ^^(u * s) - E T? B (9im (u * rj) \\ 2 ds, 

J\o,i] m 
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where r] = (771, ... , r] m ), and rji,...,r] m are independent Bernoulli random variables 77^ £ Be(sj), j = 1, . . . , m. 
Then &a, m (u) is the squared IMSE of approximation for Bp sTn (u * t),t £ [0, l] m , by the MPLI with 2 m 
observations in the vertices of unit hypercube. 

In the following theorem, we provide an exact asymptotics for the IMSE of a local stationary field approx- 
imation by MPLI when a cross regular sequence of sampling designs is used. 

Theorem 1 Let X £ 6"(X>, c(-)) 6e a random field approximated by the MPLI Xn(X,Tn), where T/v is 
cRS(h,7r,l). Then 

k 

Vi 



X - X N I ~ -£r > as N -> 00, 



1 'V 



w/iere 

Vj - = / Ci (t)6 Q .i.(i^(t^))dt>0, 
andHj(V) := {l/h j (t Lj _ 1+1 ),...,l/h j (t Lj )) > j = l,...,k. 

Remark 1 If for the j-th component, the uniform withincomponent knot distribution is used, i.e., hj(s) = 1, 
s £ [0,1], then the asymptotic constant is reduced to 

Vj = b a ,u \ Cj(t)dt, 
Jv 

where b a . t i. := 6^^.(1^). 

In Theorem [H the approximation accuracy is determined by the sampling grid sizes rij. The next theorem 
provides the asymptotically optimal inter component knot distribution for a given total number of observation 
points N. Denote by 

k \ I v-^ 1 \ _ TT h/ a 



1 / "7 



where d-p is the harmonic mean of the smoothness parameters o^, j = 1, . . . , d. 

Theorem 2 Let X £ i3f*(D,c(-)) 5e a /oca? stationary random field approximated by the MPLI X^{X,T^), 
where Tn is cRS(h,ir,l). Then 

\\X-X N \\l > k ^ asiV^oo. (5) 
Moreover, for the asymptotically optimal intercomponent knot allocation, 

n ho P t r -j- 1 as N -)• 00, j = 1, . . . , fc, (6) 

i/ie equality in © is attained asymptotically. 

The above result agrees with the intuition that more points should be distributed in directions with lower 
smoothness parameters. Note that the optimal intercomponent knot distribution leads to an increased ap- 
proximation rate. 
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Remark 2 Let X G Bf^iV, c(-)) with k = d and a.i ^ aj for some i,j = 1, . . . , d, and a := mirij = i ... d ai, i.e., 
p > a. Consider the approximation with uniform intercomponent knot distribution, n\ = ■ ■ ■ = ~ N l l d . 
Then by Theorem^ we have 

II X — Xn || 2 



jy«/(2d) • 

On the other hand, the sampling distribution (|6|) gives 

1 1 



X-X N \\ 2 x — — < 



jyp/2 i\Ta/(2d) ' 

Example 3. Let d = k = 2, a\ = 2/3, a<i = 5/3. Then for n\ = ri2, the approximation rate is 

N -a/2d = jy-1/6 while 

using the asymptotically optimal intercomponent distribution we obtain the rate 

N - P /2 = jV-1/4.2 < /^-l/6_ 

In general setting, numerical procedures can be used for finding optimal densities. However, in practice 
such methods are very computationally demanding. We present a simplification of the asymptotic constant 
expression for one-dimensional components. Further, in this case, we provide the exact formula for the 
density minimizing the asymptotic constant. For a random field X 6 B" (T> , c(-)) , define the integrated local 
stationarity functions 



Cjy 



/(*£,•):=/ Cj(t)dti . . .dt L 1 dt L+ i...dt d , t L e[0,l], j = l,...,k. 

■/[0,l]<*-i 

Moreover, for < 8 < 2, let 

2 1 
0/3 : " (B + 1)(B + 2) "6" 

Proposition 1 Let X £ ^"(X), c(-)) be a random field approximated by the MPLI Xn(X,Tn), where Tjy is 
cRS(h,n,l). If for some j, 1 < j < k, lj = 1, then for any regular density hj(-), we have 

Vj = a aj [ Cj (t L ) hj it,., )~ a i dt L . 
Jo 

The density minimizing Vj is given by 

ho P t(t L] ) = r i ^f^' . e [0, 1], 

Jo OS"/, ) </-/.. 

where jj := 1/(1 + aj). Furthermore, for such chosen density, we get 

v j,opt = 1 1 Cj I l 7j • 

In the subsequent proposition, we give an upper bound for the approximation error together with expres- 
sions for generating densities minimizing this upper bound, called suboptimal densities. 

Proposition 2 Let X £ B"(V,c(-)) be a random field approximated by the MPLI Xn(X,Tn), where Tjy is 
cRS(h,TT,l). Then 

k 

\\X -X N \\ 2 2 <Y^^i- as N ^00, 

3=1 3 
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where 



Wj = l] +aj/2 (a aj + ^ J )/,,(/, "xll. j = 1,... ,k. 



The density minimizing Wj is given by 



hj^uboptitL,) - i — L — , t Lj £ [0,1], 

Jo ( j-,. Y <h,., 

where jj := 1/(1 + a.j), j = 1, . . . , k. Furthermore, for such chosen densities, we get 

w j,subopt = lj ^ (® a i 1 1 ^3 I I -yj I J = 1) ■ ■ • ' ^ - 

Now we focus on random fields satisfying the introduced Holder type conditions. In this case, we provide 
results for the uniform mean square norm of approximation error — X^H^. The following proposi- 
tion provides an upper bound for the accuracy of MPLI for Holder classes of continuous and continuously 
differ entiable fields. 

Proposition 3 Let X G be a random field approximated by the MPLI Xn(X,Tn), where XJv is 

cRS(h,7r,l). 



i= i v 



(i) IfXe C?(V,C), then 



for positive constants c% , . . . , cp. . 
(ii) IfXe Cl' a {V,C), then 

for positive constants d%, . . . , dk . 
Remark 3 It follows from the proof of Proposition [3] that ([7]) holds if 

< : ] 2 "■/;•"-/>;■. , 1 /,. 

where Dj := 1/ min se j 0> i] hj(s), j = 1, . . . , k. Therefore the constants depend only on the parameters of the 
Holder class and the corresponding sampling design. Similar formulas can be obtained for d\, .. . ,dk in ([8]). 

In addition, we provide the inter component knot distribution leading to an increased rate of the upper bounds 
obtained in Proposition [3j 

Remark 4 Let X € C(T>) be a random field approximated by the MPLI Xn(X, Tjv), where T/v is cRS(h, ir, I). 

(i) If X € C?(V, C) and n 3 ~ iV^M , j = l,...,k, where p = (Ei=i h/ati)' 1 , then 

\\ x ~ x n\\ 00 = 0{N- po/2 ) as N -> oo. 

(ii) If X eC\' a {V,C) andrij ~ Arw/( 2 +«i),j = l,...,jfe, where Pl = (£i=i + a.)) -1 , then 

\\ x ~ x n\\ 00 = 0{N- pl/2 ) as N -> oo. 
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The approximation rates obtained in the above remark are optimal in a certain sense, i.e., the rate of con- 
verge nce can not be improved in general for random fields satisfying Holder type condition (see, e.g., iRitterl . 



200(3 ) . Moreover, these rates correspond to th e optimal approximation rates for anisotropic Nikolskii-Holder 



classes (see, e.g., 
classes. 



Yanjie and Yongpingj . l2000l ) , which are deterministic analogues of the introduced Holder 



3 Numerical Experiments 

In this section, we present some examples illustrating the obtained results. For given knot densities and co- 
variance functions, first the pointwise approximation errors are found analytically. Then numerical integration 
is used to evaluate the approximation errors on the entire unit hypercube. Let 

6 N (h, vr)(t) = 5 N (X, X N , T N (h, vr, l))(t) := X(t) - X N (X, T N {h, n, l))(t), t G [0, 

be the deviation field for the approximation of X by the MPLI with N knots, where Tjy is cRS(h, ir, 1), and 
write 

e N (h,ir) := \\ 6 N (h, 7r) 1 1 2 

for the corresponding IMSE. We write h un i(-), to denote the vector of withincomponent uniform densities. 
Analogously, by vr ura j(-) we denote the uniform interdimensional knot distribution, i.e., n\ = . . . = n^. 
Example 4. Let V = [0, l] 3 and 

X(t) = B a>l (t), 

where a = (1/2,3/2) and 1 = (1,2). Then X G Bf([0, l] 3 , c(-)), with c(t) = (1,1), t G [0,1] 3 , k = 2, ct* = 
(1/2,3/2,3/2). We compare behavior of aNihunii^uni 

) and eN(h un i,TTopt), where ir opt given by Theorem [2j 
Observe that by using the asymptotically optimal inter component distribution, we obtain a gain in the rate 
of approximation. Figure CD shows the (fitted) values of the squared IMSEs e 2 N {h un i,ir un i) and e 2 N (h un i, 7r opt ) 
in a log-log scale. In such scale, the slopes of fitted lines correspond to the rates of approximation. These 

-1.5 

-2 
-2.5 

-3 
-3.5 

-4 
-4.5 

~ 5 4 5 6 7 8 9 10 11 12 13 14 

log(N) 

Figure 1: The (fitted) plots of e 2 N (h un i,TT un i) (solid line), e 2 N {h un i,ir op t) (dash line) versus iV in a log-log scale. 

plots represent the following asymptotic behavior: 

e 2 N (h uni ,TT um ) ~ 0.3667iV" 1/6 + 0.0935iV" 1/2 ~ 0.3667iV~ 1/6 , 
e 2 N {h um ,7r opt ) ~ 0.4245iV- 3/10 as iV -> oo. 




S 



Example 5. Let V 

function 



[0, l] 2 and define X(t) 



Cov(X(t),X(s)) 



X(ti,t2) to be a zero mean Gaussian field with covariance 
1 1 , 



+ 0.1 



+ 0.1 



■cxp( 



Then X G SfQO, l] 2 , c(-)) with c(t) = ci(t) = 2/(||t|| 2 + 0.1) 2 , t G [0, l] 2 , a = 1, a* = (1,1), 1 = 2, and 
k = 1. The field has one component, hence the uniform interdimensional knot distribution is used. Theorem 
[2] provides the formula for the suboptimal withincomponent density. Figure [2ja) shows the (fitted) values 
of the squared IMSEs e 2 v (/i un j, ir un i) and e 2 N {h su b pt, vr un j). Figure [2(b) demonstrates the convergence of the 
scaled squared approximation error ^^'^G^{hsubopti^uni) to the asymptotic constant obtained in Theorem [Tl 
Note that utilizing the suboptimal withincomponent density leads to a significant reduction of the asymptotic 



6 
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5 
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4 
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3 
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4.5 5 
log(N) 



400 500 
N 



(a) 



(b) 



Figure 2: (a) The (fitted) plots of e 2 N (h 

unii^un-i) (dashed line) and eff(h su i, pt,TT un i) (solid line) versus N in 
a log-log scale, (b) The convergence of N ^e 2 N {h su b O pt^uni) (solid line) to the asymptotic constant (dashed 
line) . 

constant, as compared to the uniform withincomponent knot distribution. 



4 Proofs 



Proof of Theorem [7J First we investigate the asymptotic behavior of the approximation error ejy(t) := 
1 1 X(t) — Xjv(t) 1 1 for any t G T>\, i G I, where I := {i = . . . , id), < i\. < n£ — 1, k = 1, . . . , d}, when the 
number of knots ./V tends to infinity. Further, we find the asymptotic form of the IMSE 

e N ■= (^j^e N {tfdt^j ' 

for any positive continuous densities /ii(-), . . . , hk(-)- We start by observing that 

e^(t) 2 = E(X(t) - X N (t)) 2 = E(E r ,(X(t i + r^) - X(t))) 2 
= E v4 E ((X(t; + r^) - X(t)) (X{U + r^) - X(t))) 



1 



E^E ((X(t ; + ri *r7)-X(t)) 2 +(X(ti + ri *£)-X(t)) 2 -(X(ti + r&ny-Xfa + n**)) 2 ) , (9) 



where £ is an independent copy of rj. Further, the property ([2]) together with the uniform continuity and 
positiveness of local stationarity functions Ci(-), . . . , Cfc(-) imply that 



ejv(t) 2 



/ k 

- ( ^^(tOE^ (||rf'*(77 



+ 



r?*(£' 




l + ?jv,i(t)), (10) 
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where £n '■= max{|(7jv,i(t)|, t S A, i S 1} = o(l) as N — > oo (cf. ISeleznievl . l2000l ) . It follows from the definition 
and the mean (integral) value theorem that 



1 



1 1 \ 

" " h* d (w d , d )n*J ' W M S -^'M^ J — 1, ■ ■ ■ ,d. 



Denote by Wi := (wi^, . . . ,Wd,i d )- Now the definition of cRS(h, tt,1) implies 
/ 1 1 \ 1 



njhj(w L ._ 1+1>i )' " ' ' rijhjiwL^i ) j nj 



flj(wf), j = l,...,k, 



where Hj(t j ) := (l/hj(t Lj _ 1+ i), l/hj(t Lj )), j = 1, . . . , k. Consequently, 

/ k 



J\\\ a i 



-WHjiw!) * -£ 3 )\\ aj ) (l + o(l)) asiV^oo. 



Applying the uniform continuity of h(-) yields 

1 / k 

e N (t) 2 = - j2 n P c ^ E ^(\\ H ^i) * - sj )\\ aj + mti) * « J ' - s " 
-n^(to*(^-^-)ir) )(i+o(i)) 



J2rij Hj(ti))j (1 + o(l)) as N -> oo, 



where 

C«., {j (s^(t{)) :=iE^(||^(tO * (V - + Hfl-^tf) * (£ j - s^ll^ - H^-CtJ) * (V J ~ &)r 

= | B a j.illjW!) * s j ) - E v B aj!h (H 3 (ti)* V i) |[ . 

Let T>\ = T>\ x • • • x 2?^ and denote by \T>\\ the volume of hyperrectangle T>\. Then 

e ^ = EX tN(t) 2 dt= (Y,J^ ^n7^ Cj (tOC a ., z .(^;H i (tf))^(l + o(l)) 

= (EE^^ti) / ^(^;^i(tl))^Pi|)(l + o(l)) 

V iel 3=1 Jvj J 

\ j'=i iel / 
Now the Riemann integrability of the functions Cj(t)b a . t i.(Hj(t 3 )), j = 1, . . . , k, gives 



E n P O f ^v/.. (^i(^))*J (1 + °(!)) = y E ^7 J C 1 + °(!)) as iV ^ oo. 
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Note that for any u E , bp %m (xi) > 0, otherwise the fractional Brownian field is degenerated (cf. ISeleznjevI , 
2000). Consequently, Vj > 0, j = 1, . . . , k. This completes the proof. 



Proof of Theorem O Note that by the inequality for the arithmetic and geometric means, 

l/k 



n/ 



3=1 3 \j=l 3 



with equality if only if 



-l _ u 3 



V 



j = l,...,k. 



Hence, the equality is attained for hj = [yvj) 1 l a i , j = 1, . . . , k. Let 

rij = \hj\ ~ iyvj) 1 ' aj as N — > oo. 
The total number of observations satisfies 

d k 

N = (nj + 1) • • • (n* d + 1) ~ Yl n * = II n *j = M as N ~~> °°- 

i=l j=l 

This implies that for the asymptotically optimal intercomponent knot distribution 

it 



N ~ M ~ v 1/p Yl 

3=1 



and therefore, 



i/ ~ iV p K _p as iV oo. 
By equation (llip . the asymptotically optimal intercomponent knot distribution is 

jSfp/ocjy 1 / 01 ] 



K P/<*3 



as iV — > oo, j = 1, . . . , k. 



(11) 



Moreover, with such chosen knot distribution, the equality in ([5]) is attained asymptotically. This completes 
the proof. 



Proof of Proposition The proof is a straightforward implication of the assumptions and equation (fTUj) . The 
exact constant and the expression for the optimal density are due to lSeleznjevi ([2000 ). 



Proof of Proposition [H The first steps of the proof repeat those of Theorem [TJ By (I10p , we have 



ejv(t) 



Vi=i 

<5(E^(ti)E^(|k*(^-^) 
\j=i 



^c^tOE^ r/*(V s' 



+ rf*(£ J -s J ) - r^*(r/ J -^) 
+ ||r?*(^-^')|p') J(l + o(l)) 
(1 + o(l)) as iV -> oo. 



(l + o(l)) 
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For any nonnegative numbers a\, ■ ■ . , and any a £ M+ , the inequality 



/ k \ a k 

E«* ^ Q E 

\i=l / i=l 



a? 



(12) 



holds, and consequently, 



e N (t) 2 < J> W /2 £ E^r^l^ - Sm |)^ (1 + o(l)) 

\j=l m=Lj_i+l / 

= fE^W^ E ^%((l-^r^m + (l- S m)C))(l + o(l)). 
\j=l m=Lj_i+l / 

By the mean value theorem and the uniform continuity of withincomponent densities, we obtain 

/ k Lj \ 

e7v(t) 2 < Yl c ^ l T E ( h i(ti,m)r aj ((1 " s m )^s m + (1 - 8 m )s%) (1 + o(l)) as N -> oo. 



m=Lj- i + l 



Proceeding now to the calculation of the IMSE, we get 



(fc -kj 
£ £ C,(t,)jJ* ' V E (^(^m))^ 7- 
iel j=l m=L,_i + l ^ •? 



+ !)(«,- + 2) 



IAI (1 + 0(1)), 



where 



Jo 



(atj + l)(ay + 2) 

Now the Riemann integrability of Cj(t)hj(t m )~ aj , j = 1, . . . , k, together with the definition of integrated local 
stationarity functions imply that 



e%< (j2^j l f 2 ( a "i + l) E E c ^ t o(^(*i,m))" aj iAi)(i+o(i)) 

\i=i n j ^ ' m=ij_i+l iel / 



i{t m )- a idt (l + o(l)) 



E -Tv-/'" 2 + jf OC/. >;(.'/..; "-,//, ^ (1 + o(l)) as N -> oo. 



The expression for the suboptimal density is due to ISeleznievI (|2000l ). This completes the proof. 



Proof of Proposition [3 We start by proving (*). Let X G Cf([0, l] d ,C) and consider t € V h i € I. Applying 
the Holder condition ([1]) to equation ([9]) yields 

ejv(t) 2 = Ie^E ((X(ti + r^) - X(t)) 2 +(X(t; + ri *0 - X(t)) 2 - (X(t, + r^r?) - X(t 5 + r;*^)) 2 ) 
<CE v \\r i *r } \\ a = CE v J2\\^ j *^\\ aj <C^/; j/2 £ E,, (n, m \ Vm - s m \) a > , 

3=1 j=l m=ij_i+l 
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« 1 



where the last inequality follows from (|12p . Furthermore, since max sg [ 0j i] ((1 — + (1 — s)s aj ) =2 
we obtain 

e N (tf<cj2if 2 jr rj((i- 8m r Sm +(i- Sm >s)<^2-f/ 2 £ r.%. 

j'=l m=Lj_i+l J = l m=Lj_i+l 

By the regularity of the generating densities, we have that r^ m < \ j{n* m min s6 r 0) i] h* m {s)), i € I, m = 1, . . . , d. 
Moreover, the definition of cRS(h, tt, 1) implies the following uniform bound for the squared approximation 
accuracy 

\\X-X N \t = max4(t) < c£ 2 -^ 2 (^J 
with Dj = 1/ min se [g hj(s), j = 1, . . . ,k. Finally, we obtain the required assertion 

\\x-x»\L<JcT,-%j2> 

3=1 V 

where c] := 2 "•/!'"• 2 D°? > 0, j = 1, . . . , k. 

For the smooth case, we use the multivariate Taylor formula to obtain the following representation of the 
deviation field 

6 n (t) := X(t) - X;v(t) = E,, ^jf 5Z^(ti + nr ; * (77 - s))r ij (r ?i - Sj )du \ , teB i) t = t i + s*r i) 

where 77 = (771, . . . , 77^) and 771, . . . , r] d are independent Bernoulli random variables, 7/j G Be(sj), j = 1, . . . , d. 
Introducing an auxiliary uniform random variable U E ^(0, 1) we get 

6 n (t) = J2 E v,u {X'j(ti + U(rj - s)*r i )r iJ (77 j - s,)) 
i=i 
d 

= ^E^c/^X^i + Z7(r?i — si)ri,i, . . . ,iy + f(r?j - Sj)ry, . . . ,\ d + U(r) d - s d )r ijd ) 
i=i 

- Xj^i + J7(77i - si)^!, . . . . . . ,t ijrf + U(r] d - s d )r i:d fj (77,- - Sj), 

since for any j = 1, . . . , d, 

E v (X' j (t it i + U(rji - si)r i;1 , . . . ,iy, . . . + U(r] d - s d )r- hd )(rij - Sj)) 

= ^r n ,...,r, j - 1 , Vj+1 ,...,rid{Xj(ti tl + U('q 1 - s^r^i, . . . ,t lt j, . . . , + U(j] d - s d )r i)d )E^ (77,- - Sj)) = 0. 
The triangle inequality and the condition @ imply that 

e N (t)<^2VCV 3 r}p /2 , 
i=i 

for some positive constants Vj,j = l,...,d. Analogously to (i), the required assertion follows from the 
regularity of the generating densities and the definition of cRS(h, tt, 1). This completes the proof. 
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